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CN , ABSTRACT 

il We report the discovery of planet Kepler-12b (KOI-20), which at L695± 0.030 Rj is among the handful 

D , of planets with super-inflated radii above L65 Rj. Orbiting its slightly evolved GO host with a 4.438-day 

C/^ ■ period, this 0.431 ± 0.041 Mj planet is the least-irradiated within this largest-planet-radius group, which has 

I important implications for planetary physics. The planet's inflated radius and low mass lead to a very low 

density of 0.1 1 1 ± 0.010 g cm"^. We detect the occultation of the planet at a significance of 3.7a in the Kepler 
bandpass. This yields a geometric albedo of 0.14 ±0.04; the planetary flux is due to a combination of scattered 
Oh \ light and emitted thermal flux. We use multiple observations with Warm Spitzer to detect the occultation at 

. la and 4(t in the 3.6 and 4.5 p.m bandpasses, respectively. The occultation photometry timing is consistent 

with a circular orbit, at e < 0.01 (Icr), and e < 0.09 (3a). The occultation detections across the three bands 
favor an atmospheric model with no dayside temperature inversion. The Kepler occultation detection provides 
significant leverage, but conclusions regarding temperature structure are preliminary, given our ignorance of 
O ' opacity sources at optical wavelengths in hot Jupiter atmospheres. If Kepler-12b and HD 209458b, which 

intercept similar incident stellar fluxes, have the same heavy element masses, the interior energy source needed 
C/2 , to explain the large radius of Kepler-12b is three times larger than that of HD 209458b. This may suggest that 

^ , ■ more than one radius-inflation mechanism is at work for Kepler- 12b, or that it is less heavy-element rich than 

other transiting planets. 

Subject headings: planetary systems; stars: individual: (Kepler-12, KOI-20, KIC 1 1804465), planets and satel- 
lites: atmospheres, techniques: spectroscopic 
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1. INTRODUCTION 

Transiting planets represent an opportunity to understand 
the physics of diverse classes of planets, including mass- 
radius regimes not found in the solar system. The knowl- 
edge of the mass and radius of an object immediately yields 
the bulk density, which can be compared to models to yield 
insight into the pla net's internal co mposition, temperature, 
and structure (e.g.. Miller & Fortneyll20 11). Subsequent ob- 
servations, at the time of the planet's occultation (secondary 
eclipse) allow for the detection of light emitted or scattered 
by the planet's atmosphere, which can give clue s to a planet's 
daysi d e temperature structure and chemistry jMarlev et alj 
I2007t ISeager & DemingI 12010 ). NASA's Kepler Mission 
was launched on 7 March 2009 with the goal of finding 
Earth-sized planets i n Earth-like orbits around Sun-like stars 
(iBorucki et al.ll2010l) . While working towards this multi-year 
goal, it is also finding an interesting menagerie of larger and 
hotter planets that are aiding our understanding of planetary 
physics. 

Early on in the mission, followup radial velocity resources 
preferentially went to giant planets, for which it would be 
relatively easy to confirm their planetary nature through a 
measurement of planetary mass. This is how the confirma- 
tion of planet Kepler- 12b was made, at first glance a rel- 
atively standard "hot Jupiter" in a 4.438 day orbit. How- 
ever, upon further inspection, the mass and radius of Kepler- 
12b make it an interesting planet from the standpoint of the 
now-familiar "radius anomaly" of transit ing giant planets (e.g. 
Charbonneau et al.ll2007HBurrows etaLl l2007; La ughlin et alj 
201 Ih . Given our current understanding of strongly-irradiated 
giant planet thermal evolution, around 1/3 to 1/2 of known 
transiting planets are larger than models predict for several- 
Gyr-old planets that cool and contract under intense stellar 
irradiation (Mil ler et alJl2009l) . 

The observation that many Jupiter- and Saturn-mass planets 
are be larger than 1.0 Jupiter-radii can be readily understood. 
It is the magnitude of the effect that still needs explanation. 
The first models of strongly irradiated planets yielded the pre- 
diction that these close-in plan ets would be i nflated in radius 
compared to Jupiter and Saturn (iGuillot et al.ll996 ). The high 
incident flux drives the radiative convective boundary from 
less than a bar, as in Jupiter, to pressures near a kilobar. The 
thick radiative zone transports less flux than a fully convective 
atmosphere, thereby slowing interior cooling, which slows 
contraction. A fairly uniform prediction of these strongly irra- 
diated models is that \ .2-\3R} is about t he largest radii pre- 
dicted for planets several gigavears-old ([ Bodenheim er et al. 
2003; BuiT ows et alj|2007l: iFortnev et alj |2007: Baraff e et al. 
2008i) . However, planets commonly exceed this value. 

The mechanism that leads to the radius anomaly has not yet 
been definitively identified. However, constraints are emerg- 
ing. One is planet radius vs. incident flux, which could also 
be thought of as radius vs. equilibrium temperature, with 
an assumption regarding planetary Bond albedos. Figure [T] 
shows planet radii vs. incident flux for the transiting systems 
with confirmed masses. Since low-mass planets are relatively 
easier to in flate to large radii than higher mass planets (e.g. 
iMifler et~l . 2009), we plot the planets in three mass bins. The 
lowest mass bin is Saturn-like masses, while the middle mass 
bin is Jupiter-like masses. The upper mass bin ends at 13 
Mj, the deuterium burning limit. Kepler-12b is shown as a 
black filled cuxle. The l argest radius planets are generally 
the most highly irradiated (iKovacs et alJl2010l : lLaughlin et alJ 



l20Tlt lBatvgin"etani20Tll) . The near-universality of the in- 
flation, especially at high incident fluxes, now clearly argues 
for a m echanism that affects all close-in planets (IFortnev et alj 
l2006h . rather than one that affects only some planets. The dis- 
tribution of the radii could then be understood in terms of dif- 
fering magnitudes of the inflation mechanism, together with 
different abundances of heavy elements within the planets 
("Fortnev et al."2006!: 'Guillot e t alj|2006t iBurrows et alj|2007l: 
Miner & Fortney 201 1; Batvgi n et al.ll201 Ih . 

Within this emerging picture, outlier points are particularly 
interesting: those that are especially large, given their incident 
flux. These are the super-inflated pl anets with r adii of 1 .7 R} 
or larger. These include WA SP-12b dHebb et al.l| 2009), TrES- 
4b (Mandushevetal. 2007t ISozzetti et al.ll2009l) rWASP-17b 
(lAnderson et ah i2010 ), and now Kepler-12b, which is the 
least irradiated of the four. In the following we describe the 
discovery of Kepler- 12b, along with the initial characteriza- 
tion of the planet's atmosphere. 

Transiting planets enable the characterization of exoplanet 
atmospheres. The Spitzer Space Telescope has been espe- 
cially useful for probing the dayside temperature structure 
of close-in planetary atmospheres, as thermal emission from 
the planets can readily be detected by Spitzer at wavelengths 
longer than 3 pva. Data sets are becoming large enough that 
one c an begin to search for correlations in the c urrent detec- 
tions dKnutson et al.ll20ToHCowan & AgolllIoTlh . 

A powerful new constraint of the past two years is the possi- 
bility of joint constraints in the infrared, from Spi tzer, and the 
optical, fro m space t elescopes like CoR oT (e.g.,|G illon et ajj 
[2010 ; D eming et all 12011) and Kepler (iPesert et al . 201 Ijj). 
The leverage from optical wavelengths comes from a mea- 
surement (or upper limit) of the geometric albedo of the 
planet's atmosphere, although this is complicated by a mix 
of thermal emission and scattered light both contributing 
for these planets. Detection of relatively low geometric 
albedos Ag < 0.15 is c onsistent with cloud-free models of 
hot Ju piter atmospheres dSudarskv et alj|2003l : IBurrows et alj 
l2008h . and can inform our understanding of what causes 
the temperature inversions in many hot Jupiter atmospheres 
(Spiegel & Burrows 2010). 

In this paper we discuss all aspects of the detection, valida- 
tion, confirmation, and characterization of the planet. Section 
|2] discusses the detection of the planet by Kepler, while ^ 
covers false-positive rejection and radial velocity confirma- 
tion. Section m gives the global fit to all data sets to derive 
stellar and planetary parameters, while ^ concerns the ob- 
servational and modeling aspects of atmospheric characteri- 
zation. Section|6]is a discussion of the planet's inflated radius 
amongst its peers, while Ogives our conclusions. 

2. DISCOVERY 

The Kepler science data for t he primary transit sea rch mis- 
sion are the long cadence data dJenkins et alJl2010 F). These 
consist of sums close to 30 minutes of each pixel in the aper- 
ture containing the target star in question. These data pro- 
ceed through an analysis pipeline to produce corrected pixel 
data, then simple unweighted aperture photometry sums are 
formed to p roduce a photometric time series for each object 
dJenkins et a l. 2010c). The many thousands of photometric 
time series are then pro cessed by the transiti ng planet search 
(TPS) pipehne element dJenkins et alj|2010cl) . 

The candidate transit events identified by TPS are also vet- 
ted by visual inspection. The light curves produced by the 
photometry pipehne tend to show drifts due to an extremely 
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small, slow focus change dJenkins et al.ll20103) . and there are 
also sometimes low frequency variations in the stellar signal 
that can make analysis of the transit somewhat problematic. 
These low-frequency effects can be removed by modest fil- 
ters that have only an insignificant effect on the transit signal 
(iKoch e t al. 2010). The unfolded and folded light curves for 
Kepler- 12b produced in this manner are shown in Figure |2] 

C entroid analysis wa s performed using both difference im- 
age ([Torres et al.l201 ll) and photocenter motion (iJenkins et alj 
12010a') techniques using Ql through Q4 data. This analysis 
indicates that the object with the transiting signal is within 
0.01 pixels (0.04 arcsec) of Kepler-12, which is the 3cr ra- 
dius of confusion (including systematic biases) for these tech- 
niques. 

The parent star, Kepler Input Catalog (KIC) identification 
number 11804465, has a magnitude in the Kepler band of 
13.438. The KIC used ground-based multi-band photometry 
to assign an effective temperature and surface gravity of Teff 
= 6012 K and log g = 4.47 (cgs) to Kepler-12, corresponding 
to a late-F or early-G dwarf. Stellar gravities in this part of 
the H-R diagram are difficult to determine from photometry 
alone, and one of our conclusions based on high-resolution 
spectroscopy and light curve analyses in §4 is that the star is 
near the end of its main-sequence lifetime, with a radius that 
has expanded to R^, = 1.483 ±0.027 Rq and a surface grav- 
ity of log g = 4.175 ± 0.013. In turn, this implies an inflated 
radius for the planet candidate, originally known as Kepler 
Object of Interest (KOI)-20 (Borucki et al. 201 1). This con- 
clusion is hard to avoid, because the relatively long duration of 
the transit, more than 5 hr from first to last contact, demands 
a low density and expanded radius for the star. 

3. CONFIRMATION: FOLLOW-UP OBSERVATIONS 
3.L High Resolution Imaging from Large Telescopes 

Blends due to unresolved stellar companions (associated or 
background) can only be ruled out with direct imaging from 
large telescopes. In Figure[3]we show an image of Kepler-12 
taken with the Keck I telescope guide camera, showing 9x9 
arcsec taken in 0.8 arcsec seeing. This 1.0 second exposure 
was taken with a BG38 filter, making the passband roughly 
400 - 800 nm, similar to that of Kepler. Contours show sur- 
face brightness relative to the core. No companion is seen 
down to 7 magnitudes fainter than Kepler-12 beyond ^1 arc- 
sec from it. Thus, there is no evidence of a star that could 
be an eclipsing binary, consistent with the lack of astrometric 
displacement during transit. 

In addition, speckle observations using the WIYN tele- 
scope were made on the night of 18/19 June 2010, as part of 
the Kepler foUowup program of S. Howell and collaborators 
dHowell et al.ll201 lb . No additional source were seen to 3.69 
magnitudes fainter in /?-band and 2.17 magnitudes fainter in 
V-band in an annulus around the star spanning between 0.1- 
0.3 arcsec in radius. No companions could be seen as close as 
the diffraction limit (0.05 arcsec from the star) or as far as the 
edge of the 2.8 x 2.8 arcsecond FOV. 

Near-infrared adaptive optics imaging of Kepler-12 was ob- 
tained on the night of 08 September 2009 UT with the Palo- 
mar Hale 200in telescope and the PHARO near-infrared cam- 
era ([Havward et al. 2001) behind the Palomar adaptive optics 
system " (Trov et al. .2000) . PHARO, a 1024 x 20124 HgCdTe 
infrared array, was utilized in the 25.1 mas/pixel mode yield- 
ing a field of view of 25". Observations were performed in 
J filter (Ao = 1.25/xm). The data were collected in a stan- 



dard 5 -point quincunx dither pattern of 5" steps interlaced 
with an off-source (60" East) sky dither pattern. Data were 
taken with integration times per frame of 60 sec (15 frames) 
for a total on-source integration time of 15 minutes. The in- 
dividual frames were reduced with a custom set of IDL rou- 
tines written for the PHARO camera and were combined into 
a single final image. The adaptive optics system guided on 
the primary target itself and produced a central core width of 
FWHM = 0.11". The final coadded image at J is shown in 
Figured 

One additional source was detected at 5" SE and Ay w 8 
magnitudes fainter than the primary target, near the limit of 
the observations. No additional sources were detected at J 
within 7 "5 of the primary target. Source detection complete- 
ness was evaluated by measuring the median level and disper- 
sion within a series of annular rings, surrounding the primary 
target. Each ring has a width of 0.11" = 1 FWHM, and each 
successive ring is stepped from the previous ring by 0.1 1" = 1 
FWHM. The median flux level and the dispersion of the in- 
dividual rings were used to set the 4a sensitivity limit within 
each ring. The measured limits are in the /-band, but have 
been converted to limits in the Kepler bandpass based upon 
the typical ni Kenier-J = 1 -28 ± .52 mag for a magnitude lim- 
ited sample (iHowell et al.li201 ll) . A summary of the detection 
efficiency as a function of distance from the primary star is 
given in Figure |5] 

3.2. Radial Velocity 

To derive the planetary mass and confirm the planetary na- 
ture of the companion, observations of the reflex motion of the 
Kepler- 12b parent star were made. The line-of-sight radial ve- 
locity (RV) variati ons of the p arent star were made with the 
HIRES instrument (IVogt etalj |l994) on Keck I. Furthermore, 
a template spectrum observation was used to determine the 
stellar Teft, metallicity, and the initial log g, using the Spec- 
troscopy Made Easy (SME) tools. The log g value from spec- 
troscopy was 4.15 ± 0.05, considerably lower from the value 
in the KIC (4.47), but in good agreement with the value ob- 
tain from the Markov-Chain Monte Carlo analysis described 
in ^ The determined Teff is 5947 ± 100 K, with a distance 
estimate of ~'600 pc. ^ We note that the star is chromospher- 
ically very quiet. Our HIRES spectra cover the Ca II H&K 
lines, and we measure a chromospheric index, S=0.128 and 
log R'hk = -5.25, indicating very low magnetic activity, con- 
sistent with an old, slowly rotating star. 

All but the last four RVs were obtained during the first 
follow-up season, during the summer of 2009. The early 
Keck-HIRES spectra were taken with two compromising at- 
tributes. With a visual magnitude of V = 13.8, Kepler-12 was 
nonetheless observed with short exposure times of typically 
10-30 minutes, yielding signal-to-noise ratios near SNR=30 
per pixel for most spectra. Such low SNR taxes the Doppler 
code that was designed for much higher SNR, near 200. Thus 
the wavelength scale and the instrumental profile were poorly 
determined, increasing the RV errors by unknown amounts. 
Moreover, all observations except the last four were made 
with a slit only 2.5 arcsec tall, preventing sky subtraction, 
which is now commonly applied to HIRES observations of 
faint Kepler stars taken after September 2009. Moonlight cer- 
tainly contaminated most of these spectra, as the moon was 
usually gibbous or full, adding systematic errors to the mea- 
sured RVs. Thus the RVs given here contain some poorly 
known errors that depend on the intensity and Doppler shift 
of the solar spectra relative to that of the star in the frame of 
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the telescope. The velocities are given in Table[T] 

Based on experience with other faint stars similarly ob- 
served, we expect true errors close to 18 m s"' due to such 
effects, which are here included in quadrature. Orbital anal- 
yses should include such uncertainties in applying weights to 
the RVs, albeit not Gaussian errors. The largest RV outlier 
to our orbital analysis is the fourth RV in Table [T] and ap- 
pears at phase 0.4 in Figure |6] This measurement was made 
near morning twilight and may be more contaminated than 
the other measurements by sky spectrum. However, the mea- 
sured mass of Kepler- 12b is only modestly sensitive to these 
outliers; the mass of Kepler-12b increases by 7% when the 
largest RV outlier to a sinusoidal model is removed and the 
data are fit again. 

The phased radial velocity curve is shown in Figure |6l 
Since the orbital ephemeris from Kepler photometry was 
known a priori, observations were preferentially made at 
quadrature to allow the most robust determination of plane- 
tary mass with the fewest number of RV points. Observa- 
tions were also made at additional phases to allow an initial 
estimate of orbital eccentricity. The radial velocity observa- 
tions can be further analyzed for bisector variations, which 
are shown in Figure |6]:. No variation that is in phase with the 
planetary orbit is found, which supports the planetary nature 
of the companion. 

The radial velocities alone suggest a modest eccentricity, 
but a circular orbit certainly could not be eliminated with this 
data set. Since the long transit duration is the driver towards a 
large stellar radius, and hence a large planet radius, consider- 
able care was taken to understand if an eccentric orbit around 
a smaller parent star cou ld lead to the observed transit light 
curve (e.g. lBarnesll2007l) . As shown in Sections |4] and |5] the 
timing and duration of the occultation put more robust con- 
straints on eccentricity. 

4. DERIVATION OF STELLAR AND PLANETARY PARAMETERS 

4.1. Kepler photometry 

Our analysis is based on the Q0-Q7 data, representing 
nearly 1.5 years of data recorded in a quasi-continuous mode. 
Kepler data are in short- (SC) and long-cadence (LC) time- 
series, which are binnings per 58.84876 s and 29.4244 min, 
respectively, o f the same CCD readouts. Eight long-cadence 
jjenkins et al.l |2010b) and 16 short cadence (Gilliland et al. 
[20^10) datasets are used as part of this study, representing 
706,135 photometric datapoints and 516 effective days of ob- 
servations, out of which 464 days have also been recorded in 
short cadence. We used the raw photometry for our purposes. 

4.2. Data analysis 

For this global analysis, we used the implementation of the 
Markov Chain Monte-Carlo (MCMC) algorithm presented 
in iGillon et all (I2009i |2010). MCMC is a Bayesian infer- 
ence method based on stochastic simulations that sample 
the posterior probability distributions of adjusted parameters 
for a given model. Our MCMC implementation uses the 
Metropolis-Hasting algorithm (e.g., Carlin & Lewis 2008) to 
perform this sampling. Our nominal model is based on a star 
and a transiting planet on a Keplerian orbit about their center 
of mass. 

Our global analysis was performed using 213 lightcurves in 
total from Kepler. For the model fitting we use only the pho- 
tometry near the transit events. Windows of width 0.8 days 
(18% of the orbit) surrounding transits were used to measure 



the local out-of-transit baseline, while minimizing the com- 
putation time. In the analysis 101 SC time-series were used 
for the transit photometry. The 1-min cadence SC lightcurves 
yields excellent con s traints on the transit parameters (e.g., 
iGilliland et al.l l20l^: lip pingl 12010 ). Furthermore 112 LC 
time-series were employed for the occultation photometry. 
Input data to the MCMC also include the 16 R V datapoints 
obtained from HIRES described in Section 13.21 and the four 
Spitzer 3.6- and 4.5-iim occultation lightcurves described in 
Sectionlsn 

The MCMC had the following set of jump parameters 
that are randomly perturbed at each step of the chains: the 
planet/star area ratio, the impact parameter b' = a cos i/R^„ 
the transit duration from first to fourth contact, the time of 
inferior conjunction To (HJD), the orbital period P (assum- 
ing no transit timing variations), K' = K\/ 1 -e^P'/^, where 
K is the radial-velocity semi-amplitude, the occultation depth 
in Kepler and both Spitzer bandp asses and the tw o param- 
eters ^/ecosi^J and ^/esimo (And erson et al.ll201 ll) . A uni- 
form prior distribution is assumed for all jump parameters. 
Kepler SC data allow a precise determination of the transit 
parameters and the stellar limb-darkening (LD) coefficients. 
We therefore assumed a quadratic law and used ci = 2mi +M2 
and C2 = Ml -2m2 as jump parameters, where ui and U2 are 
the quadratic coefficients. Those linear combinations help 
in minimizing corre lations on the uncertainties of ui and U2 
jHolman et alj|2006h . 

Three Markov chains of 10^ steps each were performed 
to derive the system parameters. Their good mixing and 
convergence were asses sed using the Gelman-Rubin statistic 
(iGelman & RubinllT99l . 

At each step, the physical parameters are determined from 
the jump parameters above and the stellar mass. The tran- 
sit and radial velocity measurements together determine the 
planet orbit and allow for a geometrical measure of the mean 
density of the host star (p*). Using the MCMC chains, the 
probability distribution on was calculated and together 
with the spectroscopically measured values and uncertain- 
ties of Tgff and [Fe/H] are used to determine consistent 
stellar parameters from Yonsei-Yale stellar evolution models 
(iDemarque et al.l 120041) . The derived stellar Tgff and p^, pa- 
rameters, compared to stellar evolution tracks, are shown in 
Figure |7] The resulting normal distribution aroud the stellar 
mass (1.166±O.O52)M0 was then used as a prior distribution 
in a new MCMC analysis, allowing the physical parameters of 
the system to be derived at each step of the chains. 

4.2.1. Model and systematics 

The K epler transit and occulta tion photometry are modeled 
with the iMandel & Apoll (12002 1) model, multiplied by a sec- 
ond order polynomial accounting for stellar and instrumental 
variability. We added a quadratic function of the PSF position 
to this basel ine m odel for the Spitzer occultation lightcurves 
(see Section lSTt . 

Baseline model coefficients are determined for each 
lightcur ve with the Sing ular Value Decomposition (SVD) 
method (Pr ess et al.] 1 19921) at each step of the MCMC. Cor- 
related noise was a ccounted for following Winn et al. (2008); 
iGillon et"an (l2010h . to ensure reliable error bars on the fit- 
ted parameters. For this purpose, we computed a scaling 
factor based on the standard deviation of the binned resid- 
uals for each lightcurve with different time bins. The error 
bars are then multiplied by this scaling factor. We obtained a 
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mean scaling factor of 1.02 for all Kepler photometry, denot- 
ing a negligible contribution from correlated noise. The mean 
global Kepler photometric RMS per 30-min bin is 159 parts 
per million (ppm). 

4.3. Results 

We show in Table|2]the median values and the correspond- 
ing 68.3% probability interval of the posterior distribution 
function (PDF) for each parameter obtained from the MCMC. 
We present in Figure [8] the phase-folded transit photometry. 
We determine a planetary radius of 1.695;!;Q g32 ^ mass 

of 0.431;!{5 o4o Mj that produces a very low mean planetary 
density of 0. 11 l![J:[|],', g cm'l 

We measure occultation depths of 0.099 ± 0.028% and 
0. 1 16 ± 0.034% in Spitzer IRAC 3.6 and 4.5/im channels re- 
spectively, consistent at the Icr level with the specific anal- 
ysis present in Section 15.11 The LD quadratic coefficients 
derived from the MCMC are mi = 0.375 ± 0.004 and U2 = 
0.250 ±0.008. Those are in good a greement with the theoret- 
ical coefficients obtained from the lClaret & BloemenI (1201 ll) 
tables of ui = 0.366 and U2 = 0.275. 

We finally determine an occultation depth of 31 ±8 ppm 
in the Kepler bandpass, which corresponds to a geometric 
albedo Ag = 0.14 ±0.04. The geometric albedo is wavelength- 
dependent and measures the ratio of the planet flux at zero 
phase angle to the flux from a Lambert sphere at the same 
distan ce and the same cross-sectional a rea as the planet (see, 
e.g.. ,Marley etan[T999l:ISudarsky et al.ll200Q) : 




where j- is the occultation depth, a the orbital semi-major 
axis and Rp the planetary radius. 

The corresponding phase-folded occultation lightcurve is 
shown in Figure|9] The combination of Spitzer and Kepler oc- 
cultations leads to a Icr orbital eccentricity signal of e < 0.01, 
while the 3<j limit is e < 0.09. We show esinuj vs. ecoso; 
from successful MCMC trials in Figure [TO] The small al- 
lowed eccentricity removes most solutions that allow fits to 
the long transit duration with smaller stellar (and planetary) 
radii. There are two paths towards a more robust constraint 
on e. One would come from many additional RV points. An 
easier path would be additional quarters of Kepler data, which 
would yield a better determination of the occultation duration, 
which constrains esinw. All system parameters are collected 
in Table |2] 

5. ATMOSPHERIC CHARACTERIZATION AT SECONDARY ECLIPSE 

As part of Spitzer program #60028 (D. Charbonneau, PI) 
a number of Kepler-detected giant planets were observed in 
order to characterize the planets' thermal emission at 3.6 and 
4.5 fim during the Warm Spitzer extended mission. The in- 
herent faintness of the planetary targets mean some stars must 
be observed more than once for adequate signal-to-noise to 
enable meaningful atmospheric characterization. 

In addition to the measurement of the depth of the occul- 
tation (or secondary eclipse), which yields a measurement of 
the planetary brightness temperature, the timing and duration 
of the occultation constrains e, as described above. The tim- 
ing of the transit constrains e cos uj where lu is the longitude 
of periapse. The duration of the transit constrains esinui. The 
former is generally easier to measure accurately than the lat- 
ter. 



5.1. Warm Spitzer Detections 

Kepler- 12 was observed during four occultations between 
August 2010 and Januar y 2011 w ith Warm-SpitzerflRAC 
C Werner et al. 2004; Fazio et al.l2004 ) at 3.6 and 4.5 ^m. Two 
occultations were gathered per bandpass and each visit lasted 
approximately 11 h. The data were obtained in full-frame 
mode (256 x 256 pixels) with an exposure time of 30.0 s per 
image which yielded 1321 images per visit. The set of obser- 
vations are shown in Table [3] 

The method we used to produc e photometri c time s eries 
from the images is described in Desert et"al] (1201 lal) . It 
consists of finding the centroid position of the stellar point 
spread function (PSF) and performing aperture photometry 
using a circular aperture on individual exposures. The images 
used are the Basic Calibrated Data (BCD) delivered by the 
Spitzer archive. These files are corrected for dark current, flat- 
fielding, detector non-linearity and converted into flux units. 
We convert the pixel intensities to electrons using the infor- 
mation on detector gain and exposure time provided in the 
FITS headers. This facilitates the evaluation of the photo- 
metric errors. We extract the UTC-based Julian date for each 
image from the FITS header and correct to mid-exposure. We 
then correct for transient pixels in each individual image us- 
ing a 20-point sliding median filter of the pixel intensity ver- 
sus time. To do so, we compare each pixel's intensity to the 
median of the 10 preceding and 10 following exposures at 
the same pixel position and we replace outliers greater than 
3 a with its median value. The fraction of pixels we cor- 
rect varies between 0.15% and 0.22% depending on the visit. 
The centroid position of the stellar PSF is determined using 
DAOPHOT-type Photometry Procedures, GCNTRD, from the 
IDL Astronomy LibrarjQ We use the APER routine to per- 
form aperture photometry with a circular aperture of variable 
radius, using radii of 1.5 to 8 pixels, in 0.5 steps. The prop- 
agated uncertainties are derived as a function of the aperture 
radius; we adopt the one which provides the smallest errors. 
We find that the transit depths and errors vary only weakly 
with the aperture radius for all the light-curves analyzed in 
this project. The optimal apertures are found to have radii of 
2.5 pixels. 

We estimate the background by fitting a Gaussian to the 
central region of the histogram of counts from the full array. 
The center of the Gaussian fit is adopted as the residual back- 
ground intensity. As akeady seen in previous Warm-Spitzer 
observations (Deming et al. 201 1; Beerer et al. 201 1), we find 
that the background varies by 20% between three distinct lev- 
els from image to image, and displays a ramp-like behavior 
as function of time. The contribution of the background to 
the total flux from the stars is low for both observations, from 
0.07% to 1.2% depending on the images. Therefore, photo- 
metric errors are not dominated by fluctuations in the back- 
ground. We used a sliding median filter to select and trim 
outliers in flux and positions greater than 4 a. This process 
removes between 0.9% and 2.8% of the data, depending on 
the visit. We also discarded the first half-hour of observations, 
which are affected by a significant telescope jitter before sta- 
bilization. The final number of photometric measurements 
used are presented in Table |3] The raw time series are pre- 
sented in the top panels of Figure [TT] 

We find that the point-to-point scatter in the photometry 
gives a typical signal-to-noise ratio of 260 and 200 per im- 

^^ |http : //idlastro . gsf c ■ nasa ■ gov/homepage ■ html| 



6 



Fortney, et al. 



age at 3.6 and 4.5 fim respectively. These correspond to 85% 
of the theoretical signal-to-noise. Therefore, the noise is dom- 
inated by Poisson photon noise. We used a transit light curve 
model multiplied by instrumental decorrelation functions to 
measure the occultation parameters a nd their uncertainties 
from the Spitzer data as described in Desert et al.l (f201 lb). We 
compute the transit light curves with the IDL transit routine 
OCCULTSMALL from .Mandel & Agoll (120021) . In the present 
case, this function depends on one parameter: the occulta- 
tion depth d. The planet-to-star radius ratio Rp/R-^, the orbital 
semi-major axis to stellar radius ratio (system scale) a/7?*, 
the mid-occultation time Tc and the impact parameter b are 
set fixed to the values derived from the Kepler lightcurves. 

The SpitzerflRAC photometry is known to be systemati- 
cally affect ed by t he so-called pixel-phase effect (see e.g., 
ICharbonne au et af] l2005: Knutson et al. 2008). This effect is 
seen as oscillations in the measured fluxes with a period of 
approximately 70 min (period of the telescope pointing jit- 
ter) and an amplitude of approximately 2% peak-to-peak. We 
decorrelated our signal in each channel using a linear function 
of time for the baseline (two parameters) and a quadratic func- 
tion of the PSF position (four parameters) to correct the data 
for each channel. We performed^a simultaneous Levenberg- 
Marquardt least-squares fit dMarkw ardt 2009) to the data to 
determine the occultation depth and instrumental model pa- 
rameters (7 in total). The errors on each photometric point 
were assumed to be identical, and were set to the rms of the 
residuals of the initial best-fit model. To obt ain an esti mate of 
the correlated and systematic errors (P ont et al.l l2006) m our 
measurements, we use the residual pe rmutation bootstrap , or 
"Prayer Bead", method as described in lDesert et alj (120091) . In 
this method, the residuals of the initial fit are shifted system- 
atically and sequentially by one frame, and then added to the 
transit light curve model before fitting again. We allow asym- 
metric error bars spanning 34% of the points above and below 
the median of the distributions to derive the Ig uncert ainties 
for each parameter as described inl Desert et alj (1201 lal) . 

We measure the occultation depths in each bandpass and for 
each individual visit. The values we measure for the depths 
are all in agreement at the la level. Furthermore the weighted 
mean averages per bandpass of the transit depths are consis- 
tent with the depths derived by the global Monte-Carlo anal- 
ysis. 

5.2. Joint Constraints on the Atmosphere 

To model the planet's atmosphere we use a one-dimensional 
plane-parallel atmosphere code that has been widely used for 
solar system planets, exoplanets, and brown dwarfs over the 
past two decades. The optical and therm al infrared radiativ e 
transfer solvers are described in detail in iToon et al] (I1989I). 
Past applications of the model include Titan ("McKav et ajj 
[l989), Uranus (Marlev & McKav 1999), gas giant exoplanets 



(iFortnev et alJ I2006t IFortnev & Marlev 2007; Fortnev 



and brown dwarfs (Marlev et al. 1996; Burrows ( 
; iMarlev et alj|2002l; Saumon & Marlev; 120081). We use 



et aD 
etall 



the co rrelated-k method for opacity tabulation (' Goody et alj 
19891). Our extensiv e opacity database is described in 
Freedman et al.l (1200 8') . We make use of tabulations of 



chemical mixing ratios from equilibrium chemistry calcu- 
lations of K. Lqdders and collaborators (Lodders 199^ 
iLodders & Feglevll2002ll2()06l) . We use the protosolar abun- 
dances of iLoddersI (|2003[k Si nce the first detect i on of 
thermal flu x from hot Jupiters (I Charbonneau et al] 120051 ; 
iDeming eFal.. .2005.) we have used the code extensively to 



model strongly irradiated planet atmospheres and have com- 
pared model spectra to observatio ns (e.g . Fortnev et al. 2005|; 
iKnutson et al.ll2009t:lDeming et alj|201 ItiDesert et ahi201 lal) . 

Planet Kepler- 12b intercepts an incident flux of 1.1 x 10^ 
erg s"' cm"^, a value just larger than the suggeste pM /pL class 
incident flux boundary proposed by IFortnev et al.l (2008). It 
was suggested that planets warmer than this boundary (pM) 
would harbor dayside temperature inversions, while those 
cooler than this boundary would not have inversions. It is 
therefore important to understand the temperature structure 
of the planet. For Kepler- 12b we show three models in Figure 
[T2l for which we plot the planet-to-star flux ratio and dayside 
P-T profiles. In red and blue are "dayside average" models 
with incident flux redistributed over the dayside only. In green 
is a model where the incident flux is cut in half, to simulate 
efficient redistribution o f energy to the night side (see, e.g. 
IFortnev & Marlevil2007l) . The model in red has a temperature 
inversion due to absorption of incident flux b y TiO and VO 
vapor (e.g. lHubenv et al.ll2003l: [Fortnev et alT2 008), while the 
blue and green models lack inversions, as TiO/VO vapor is 
removed from the opacities. The Kepler occultation depth is 
shown at 0.65 /im (diamond), while the Spitzer detections are 
shown as diamonds at 3.6 and 4.5 yum. Model band-averages 
at these wavelengths are shown as solid circles. 

The relatively flat ratio of the 3.6/4.5 diamond poin ts gen- 
erally points to a very weak or no inversion (Knutso n et alJ 
'2010'). Looking to the optical, the green model is dramat- 
ically too dim, while the blue model nearly reaches the la 
error bar Looking at the infrared, the blue point is at the 
la 4.5 /im error bar as well. The inverted model (red) has 
approximately the same T^ff (^ 1700 K) as the blue model, 
but higher fluxes in the mid infrared and lower fluxes in the 
near-infrared and optical. The Spitzer data alone do not give 
us strong leverage on the temperature structure. Any cooler 
model with an inversion (not plotted) would yield a better fit 
to Spitzer and a worse fit to Kepler. Within the selection of 
models, the brightness of the Kepler point argues for the no- 
inversion model. The flux in the Kepler band from the blue 
model is 60% scattered light, 40% thermal emission. 

Our tentative conclusion is that the blue (no inversion, 
inefficient temperature homogenization onto the night side) 
model is preferred. However, given our ignorance of the opti- 
cal opacity in these atmospheres, this conclusion is tentative. 
The relatively deep occultation in the Kepler band argues for 
an additional contribution at optical wavelengths that is not 
captured in the model. One possibi lity is that stellar flux has 
photoionized Na and K gasses ( Fortney et alj2003) . which are 
thought to be strong absorbers of stellar light (and therefore 
diminish scattered light) in hot Jupiter atmospheres. Another 
possibility is a population of small grai ns, such as silicates , 
which could scatter some stellar flux (iMarlev et alJ [19991; 
ISeager et aL 2000; Sud arsky et all I2000I). Such clouds are 
prom inent in L-dwarf atmospheres (e.g. [Ackerman & Marlevl 
[200l . 

6. DISCUSSION 

A great number of explanations have been put forward to 
explain the inflated radii of the close-in giant planets. They 
generally fa ll into s everal broad classes, and are recently 
reviewed in iFortney^ & Nettelmann (2010[) and iBaraffe et alJ 
( [20T0I) . Some argue for a delayed contractio n, due to slow ed 
energy transport in the atmosphere (Burrow s et al.l [20071) or 
the deep interior ( Chabr ier & Baraffa |2007|). O thers sug- 
gest a variety of atmospheric affects dShowma'n & GuillotI 
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200llGuillot & Showrnanl llOOllBatvgin & Stevensonll2010t 
Arras & SocratesI 120 lOt lYoudin & Mitchell 201(| that lead 
to energy dissipation into the interior. Still others suggest 
tidal dissipation in the interior due to eccentri c ity damping 
jBodenheimer et al. 2000; Ja ckson et all l2008t iMiller et alJ 
[2009; Ibgui j&Burrows 2009). 

For Kep ler- 12b we do no t find evidence for transit timing 
variations (iFord et aT]|201 l l). The RMS scatter of transit times 
about a linear ephemeris is less than one minute and only 17% 
larger than the average of the formal timing uncertainties. 
This rules out the presence of massive non-transiting planets 
very near by or in the outer 1:2 mean motion resonance. In 
principle, a more distant non-resonant planet is possible, but 
hot Ju piters rarely have a second massive planet c lose to the 
star (iWright et al.ll2009l lIoTTt iLatham et aLllIoTl . Thus, it 
is very unlikely that the inflated radius is due to eccentricity 
damping. 

Clarity on a radius-inflation mechanism has not been 
achieved, but Figure [T] appears to argue for an explana- 
tion based on the planet temperature or irradiation level of 
the atmosphere (rather than merely on orbital separation) , 
as has been shown by other authors (iKovacs et alJ 120101 ; 
Eaughlin et al. 2011; Batygin et al. 201 1). 

If the inflation mechanism can be thought of as an energy 
source that is added to the planet's deep convective interior, 
we can readily compare the energy input needed to sustain 
the radius of Kepler-12b, compared to other planets. This is 
actually more physically motivated than the more commonly 
discussed "radius anomaly," since the power needed to inflate 
the radius by a given amount, AR, is a very s t rong f unction 
of mass. In particular. Figure 6 in Mill eFetaLl (l2009h allows 
for a comparison of input power as a function of planet mass, 
for 4.5 Gyr-old model planets with 10 cores at 0.05 AU 
from the Sun. For instance, inflating a 0.2 Mj planet by 0.2 
Rj over its expected radius value takes 1 x 10^^ erg s"', while 
for a 2 Mj planet it is 2 x 10^^ erg a factor of 2000 dif- 
ference in power for a factor of 10 in mass. This is the rea- 
son why Batygin et al. (201 1) can easily expand Saturn-mass 
planets to the point of disruption via Ohmic dissipation — a 
small amount of energy goes a long way towards inflating the 
radii of low-mass planets. 

In understanding the stru cture of Kepler- 12 b, we can 
use the mode ls described in iMiller et al.l (l2009h . which are 
adapted from 'Fortne y et al ] (120071) . In particular. Table 1 in 
[Miller et al. (2009) includes the current internal power nec- 
essary to explain the radius of several inflated planets. Plan- 
ets HD 209458b ("Henrv et al.l2000t JCharbonneau et al.l2()00l) 
and TrES-4b (Mandushev et al. 2007|) are interesting points of 
comparison. HD 209458b and Kepler-12b have similar inci- 
dent fluxes, while TrES-4b and Kepler-12b have similar in- 
flated radii. 

Since Kepler-12b and HD 209458b have comparable in- 
cident stellar fluxes (that of Kepler-12b is 14% larger), one 
could easily assume that they have sim ilar interior energy 
sources (e.g. iGuillot & Showmanll200l . For HD 209458b, 
with core masses of 0, 10, and 30 M®, incident powers 
of 1.5 X l O^'', 3.8 X 102^ and 1.6 x 10^^ erg s-\ are re- 
quired (M iller et al.l l2009). Using the planetary parameters of 
Kepler- 12b with cores of 0, 10, and 30 M^, the required pow- 
ers are substantially larger. The enhancement is generally a 
factor of three larger, with values of 4.4 x 10^^, 1 . 1 x 10^^, and 
4.2 X 10^^ erg s"', respectively. This could point to more than 
one radius inflation mechanism being at play in this planet, as 



has recently been stro ngly suggested for the massive transit- 
ing planet CoRoT-2b ( IGuillot & Havel201 lb . 

The difference between Kepler-12b and HD 209458b can 
be remedied, however, if the planets have different heavy ele- 
ment masses. In particular, both planets would require power 
levels of ~ 1.6 X 10^^ erg s"' if Kepler-12b possesses only 
~ 15 of heavy elements, while HD 209458b possesses 
30 Mq. The Kepler-12b parent star has an [Fe/H]=H-0.07, 
while for HD 209458 it is H-0.02. As recently shown by 
IMiller & FortnevI (1201 1) for the colder non-inflated planets, 
for parent stars with similar stellar metallicities, a spread from 
10-30 Mq) is reasonable. Therefore it appears that the wide 
disparity in radii between these two well-studied planets could 
alternatively be due to the differences in interior heavy ele- 
ment masses. Large diversities in heavy element abundances 
are clearly needed to explain plots like Figure[T3] where plan- 
ets of similar masses can have dramatically different radii. 

For comparison, TrES-4b at 0.93 Mj and 1.78 Rj is nearly 
twice as massive as Kepler-12b, but intercepts 2.1 times 
higher incident flux. The inflation powers at 0, 10, and 30 
range from 1.0 to 3.4 x 10^**, 8-20 times larger than 
for Kepler-12b, at the same heavy element masses. Clearly 
the required energy difference between the two models does 
not scale simply with the incident flux. As discussed in 
IMiller & Fortneyl (1201 ll) as the population of cool (req< 1000 
K) non-inflated planets grows, the heavy element mass of ex- 
trasolar gas giants can become better understood as a function 
of planet mass and stellar metallicity, which will allow for 
more robust constraints on the heavy element masses of the 
inflated planets. This will in turn allow for better estimates of 
the magnitude of the additional interior energy source within 
these planets. While Kepler-12b does not quite fit the general 
trend that the highest irradiation planets are the largest, this 
trend argues for an explanation that scales with atmospheric 
temperature. A more detailed computational understanding 
of how the visible atmosphere, deep atmosphere, and convec- 
tive interior interact and feedback on each other is now clearly 
needed. 

7. CONCLUSIONS 

We report the discovery of planet Kepler- 12b from transit 
observations by Kepler. The planet has an unusually inflated 
radius and low bulk density. At its incident flux level, the large 
radius of the planet makes it somewhat of an outlier compared 
to the general empirical trend that the most inflated planets in- 
tercept the highest incident fluxes. This may require the planet 
to have an usually low mass fraction of heavy elements within 
its interior, or that more than one radius-inflation mechanism 
is at work in its interior 

The atmosphere of the planet is probed via detections of 
the occultations in the Kepler and Warm Spitzer bandpasses. 
Given the faintness of the parent star, characterization was dif- 
ficult, but all detections were made at a level of at least 3.5cr. 
A model comparison to the data yields a best-fit model that 
lacks a dayside temperature inversion, given the relatively flat 
3.6/4.5 fim ratio of the planet-to-star flux ratios, along with 
the relatively large occultation depth in the Kepler band. Ad- 
ditional Kepler data will yield more robust constraints on the 
planet's geometric albedo, orbital eccentricity, and perhaps 
phase curve information. 
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Fig. 1 . — Observed planetary radius as a function of total incident flux for the known transiting gas giants. Planets are plotted in three colors for three different 
mass cuts. Kepler- 12b is shown as a black filled circle. Note the general trend towards smaller radii with decreased insolation. There is a probable break in slope 
at incident fluxes of 1 - 2 ~ 10* erg cm"^ s"' . Planets are taken from the compilation at: http ://www.inscience.ch/transits/ 1 



TABLE 1 

Relative radial velocity and bisector span variation measurements of Kepler- 12. 



BJD 


RV 


0"RV 


BS 


TBS 




(m s-i) 


(ms-') 


(m s-') 


(m s-i) 


2455014.91234 


-12.5 


19.3 


38.3 


6.1 


2455016.79104 


45.9 


19.3 


46.2 


6.7 


2455017.10568 
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19.1 


21.1 


3.4 


2455019.11155 


36.7 


19.4 


22.3 


4.8 


2455027.05631 


-57.9 


19.4 


21.5 


4.6 


2455041.99855 


-5.1 


19.7 


17.6 


5.9 


2455042.79449 
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19.6 


32.3 
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2455073.83165 
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19.0 


-81.8 
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68.3 


6.1 


21.0 


3.5 


2455439.75382 


-29.6 


6.0 


20.2 


3.7 


2455759.86617 


-36.6 


6.0 


31.4 


3.8 


2455761.82490 


48.5 


6.4 


27.1 


4.1 
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Fig. 2. — Time series and folded transit light curve for Kepler-12b. 
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Fig. 3. — Image of Kepler-12 taken with the Keck 1 telescope guide camera, showing 9x9 arcsec taken in 0.8 arcsec seeing. Noilh is up and east is to the 
left and the pixels are 0.30 arcsec in size. The exposure time was 1.0 sec. The detector is a Photometries CCD and the filter is a BG38, making the passband 
roughly 400 - 800 nm, similar to that of Kepler. Contours show surface brightness relative to the core. No companion is seen down to 7 magnitudes fainter than 
Kepler-12 beyond ~1 arcsec from it. There is no evidence of a star that could be an eclipsing binary. 



12 



Fortney, et al. 




+6 +4 +2 0.0 -a -4 -6 +1.0 +0.5 0.0 -0.5 -1.0 

Aa [arcsec] Aa [arcsec] 



Fig. 4. — / Palomar adaptive optics image of Kepler-12. The left image displays a 15" X 15" field of view centered on the primary target. The light image 
displays a 2" X 2" field of view centered on the primary target. The four-point pattern surrounding the central point spread function core is part of the adaptive 
optics point spread function. 




Fig. 5. — Palomar detection Umits as a function of radial distance from the primary target, Kepler-12. The filled circles represent the J-band limits and each 
point represents a step in FWHM away from the primary target centroid peak. The dashed line underneath represents the J-band limits converted to Kepler 
magnitude limits if a star were to have a nominal micgpig,-—J color. (For a magnitude limited sample, the median mf^gpig,-—J = 1.28 it 0.52). 
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Fig. 6. — (a) Orbital solution for Kepler- 12b. The observed radial velocities obtained with HIRES on the Keck Telescope are plotted together with the velocity 
curve for a circular orbit with the period and time of transit fixed by the photometric ephemeris. The radial velocities have an arbitrary zero point, (b) Velocity 
residuals from the orbital solution. The rms of the velocity residuals is 24 ms"' . (c) Variation in the bisector spans for HIRES spectra. The zero point is arbitrary 
and the RMS is 59 ms"'. 
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intervals on the stellar T^ff and p^, as determined by spectroscopy and transit model fits. The 'hook' in the evolution tracks for more massive stars with convective 
cores produces a non-uniform distribution of masses with the uncertainty boxes and can produce a degeneracy in the solution for the best-fit stellar parameters. 
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Fig. 8. — Top: Kepler-12 b phase-folded transit lightcurve with best-fit model superimposed. The data are binned in 15 minute intervals. Error bars are smaller 
than the plotted datapoints. Bottom: residuals are displayed in parts per million (ppm) scale. 
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Fig. 9. — Top: Kepler-12 b phase-folded occultation lightcurve with best-fit model superimposed. The Kepler data are binned in 15 minute increments. Bottom: 
residuals are displayed in parts per million (ppm) scale. 
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Fig. 10. — Density function of the two-dimensional esintj/ecosci; successful MCMC trials (density increases from white to black). Note the different scales 
for X and y axes. The 1, 2 and 3(t confidence domains are superimposed. The weak constraint from the occultation duration leads to a significant scatter in esincj. 
As e cos w is well constrained from the occultation timing obtained by the Kepler and Spitzer photometry, this yields orbital eccentricities as high as ~ 0.09 (3<7 
upper limit) provided the argument of periastron is close to 90 or 270 degrees. In this case, the eccentricity vector points toward or away of the observer, allowing 
a wide range of es'muj values while leaving ecosw almost unchanged. 
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Fig. 1 1. — Spitzer occultadon light-curves of Kepler-12b observed in the IRAC band-pass at 3.6 (top) and 4.5 (bottom). Top panels: raw and unbinned 
light-curves. The red solid lines cor respo nd to the best fit models which include the time and position instrumental decorrelations as well as the model for the 
planetaiy transit (see details in Sect. l5.lt . Bottom panels: connected and normalized occultation light-curve with the best fit model (in red). The data are binned 
in 25 minutes intervals (50 points). 
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Fig. 12. — Main Panel: Planet-to-star flux ratios observed by Kepler and Spitzer, shown in gray. The flux ratios from thi'ee planetary models are shown for 
comparison. There is preference towards models that have no temperature inversion (blue and green). Model ratios integrated over the appropriate bandpasses 
are shown as filled circles. The Kepler occupation point strongly ai'gues for inefficient redistribution of flux, or an additional scattering component at optical 
wavelengths. Inset Panel: Atmospheric pressure-temperature profiles for the three models. 
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Fig. 13. — M ass vs. radius for planets with "well-defined parameters," as taken from http://www.inscience.ch/transits/, but also including the Kepler- 1 1 system 
ILissauer et al . 201 1). Kepler-12b is shown as an open circle, the 2nd-lowest density planet discovered. Models (solid black curves) are taken from .Fortnev et all 
120071) . The two upper curves are for pure H-He planets, at 4.5 Gyr, at 0.02 and 10 AU from the Sun. 
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TABLE 2 
Kepler- 12 system parameters 



Parameters 



Value 



Jump parameters 



Planet/star area ratio (Rp/Rsf' 
b' = a cos i/R^ [R^,] 
Transit width [d] 
7b - 2450000 [HJD] 
Orbital period P [d] 
RV^:' [ms"' d'/3] 

Y^COSLJ 

Y^sinw 
ci = 2w[ +U2 
C2 = u\ —2u2 
Occultation deptli 

Deduced stellar parameters 



013765+<'"°'"'2'* 
174+0.011 

n 1 o<;7-!+0.00009 

5004.00835^*'''""'2 



4 4^70637+0 0000012 
t.t J / 7U J / -0.0000002 

79.2!7;i 

_n (¥11+0 054 

Oooi^oDST 
1 009+!^:^^^ 

n I yi+O.lllfi 

-U.182_„,„4 

00003 1+^' ™'"" 



Ml 

Density p^, [pq] 
Surface gravity logg* 
Mass Mi, [Mq] 
Radius R* [Rq] 
Age [Gyr] 



[cgs] 



0.274;^ 
0.354^! 
4.175! 

1-1664p54 
r /ioo+o.()25 

1 •tO-3_o.o29 
^•"-0.4 



Observed stellar parameters 



Trff 

[Fe/H] 
Vsin/ 



5947 ± 100 
0.07 ± 0.04 
0.8 ±0.5 km 



Deduced planet parameters 



RVi:[ms-'] 
btransit [R*] 
^occultation 

Toccultation " 2450000 [HJD] 
Orbital semi-major axis a [AU] 
Orbital inclination i [deg] 
Orbital eccentricity e 
Argument of periastron a; [deg] 
Density pp [g cm"'] 
Surface gravity log gp [cgs] 
Mass Mp [M,„p] 
Radius Rp [Rup] 



48.2:|-4 
Q J74+0.011 

i74;8;81i 
"•^ '^-o.ou 
5010.666!<;!S^ 
0556-*<'«''" 

00 nfMi)% 
°°-'"-0.08 

< 0.01 (Ict), < 0.09 (3(t) 
182+^^ 

111+0.011 

' -0.009 

2 57+0.04 

n 437+Stl41 

-Q.Q4Q 

1 f:Qc+0.tH8 
i •"='-'-0.032 



TABLE 3 
Warm-Spltzer observations. 



Visit AOR Wavelength Obs. Date (UT) Select, points Depth (%) Weighted. Avg. depth Tbright 



1 40251392 3.6 2010-09-06 1233 O.Ml^oo^b 

3 40250880 3.6 2010-12-26 1151 0.130!o|° 0.137± 0.020% 1597±160K 

2 40251136 4.5 2010-09-15 1160 0.108^°™ 

4 40250624 4.5 2011-01-08 1212 0.12931? 0.116±0.031% 1429±190K 



